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POLYNOMIAL INTERPOLATION AND APPROXIMATION 

IN 

T. BLOOM*, L. P. BOS, J.-P. CALVI AND N. LEVENBERG 

Abstract. We update the state of the subject approximately 20 years 
after the pubUcation of [8] . This report is mostly a survey, with a sprin- 
■ kling of assorted new results throughout. 

> 
O 

^ ■ 1. Introduction. 

OO _ 

, Let zq, ...Zn be n + 1 distinct points in the plane and let / be a func- 

tion which is defined at these points. The polynomials lj{z) = Y[k=/=j{^ ~ 
^ ! ^k) / Ylk^j{zj — Zk)i j = 0, ...,n, are polynomials of degree n with lj{zk) = 
^ ' Sj^k which we call the fundamental Lagrange interpolating polynomials, or 
^ . flip's, associated to ZQ,...,Zn- The polynomial p{z) = J^^^o f{zj)lj{z) is 
d . then the unique polynomial of degree at most n satisfying p{zj) = f{zj), j = 
0, ...,n; we call it the Lagrange interpolating polynomial, or LIP, associated 
to /, zo, •••5 ^n- If r is a rectifiable Jordan curve such that the points zq, Zn 
are inside F, and / is holomorphic inside and on F, we can estimate the er- 
ror in our approximation of / by p at points inside F using the Hermite 
Remainder Formula: for any z inside F, 

where u{z) = YYk=oi^ ~ ^k)- This elementary yet fundamental formula is 
the key to proving many important results on polynomial approximation 
and interpolation. We first recall the following result of Walsh which gives 
^ ' a quantitative version of the classical Runge theorem. 

H : 

Theorem 1. (Walsh) Let K be a compact subset of the plane such that 
C\K is connected and has a Green function qk- Let R > 1, and define 

D R := {z e C : gK{z)< log R}. 

For f continuous on K , let 

(1.2) dn{f, K) := mi{\\f — pn\\K '■ Pn polynomial of degree < n} 
where ||/ -p„||i^ = sup^g^ \f{z) -pn{z)\. Then 

limsupd„(/,ir)i/" < 1/R 

n— >oo 

if and only if f is the restriction to K of a function holomorphic in D^. 



> 
OO 
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Here, C\K has a Green function means that qk is continuous and 
subharmonic in C, harmonic in C\K with gxi^) — log 1^1 bounded as \z\ — >■ 
oo, and qk = on K. This final condition says that i^" is a regular compact 
set. 

Consider the following situation. Let {znj}-, j = 0,...,n; n = 1,2,... 
be an array of points. For each / defined in a neighborhood of this array, 
we can form the sequence of LIP's {Lnf} associated to /. Let ujn{z) : = 



— Znj). An easy consequence of Theorem [T] and (II. ip is the following. 



Theorem 2. Let K G C be compact and regular with C\K connected. Let 
{znj} be an array of points in K. Then for any f which is holomorphic in 
a neighborhood of K, we have Lnf ^ f on K (uniform convergence) if and 
only if 



lim lur, 



Z) r^ + l 



6{K) 



o9k{z) 



uniformly on compact subsets of<C\K. 



Here 6{K) is the transfinite diameter of K. In [8j, several conditions on 
the array {^nj} were discussed which imply for any / which is holomorphic 
in a neighborhood of K, we have L„/ =^ f on K. The results in this 
univariate setting are well understood. In C^, d > 1, knowledge of Lagrange 
interpolation is less complete. Let Vn denote the complex vector space of 
holomorphic polynomials of degree at most n and let 

/ n + d^ 

N{n) := dimP„ 



n 



Thus 



where {cj 



Vn = span{ei, ...,6^} 
z"'^^^} are the standard basis monomials. We let 

dnN 



N 



d + 1 



For points Ci, 
of order n as 

(1-3) 



, Ctv € 'C'^, define a (generalized) Vandermonde determinant 



VDM{Cu-XN) = det[e,{Qh 
ei(C2) .. 



,N 



det 



ei(Ci; 
. e7v(Ci 



Given N points An = {A 



e7v(C2) 
AnN} with 
VDM{Ani,...,A 



nl ) 



nN 



we can form the FLIP's 

VDM{Anl,...,X,...,AnN) 



Inj {x) 



ei(Cw) 
e7v(Civ) 

7^0, 

J 



N. 



VDMiAnl,...,AnN) 

In the one (complex) variable case, we get cancellation in this ratio so that 
the formulas for the FLIP's simplify. In general, we still have lnj{A 
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and Inj € Vn since Inj is a linear combination of ei, .., cat. For / defined at 
the points in A^, 

N 

(1.4) (L„/)(x):=^/(A„,)/„,(x) 

i=i 

is the Lagrange interpolating polynomial (LIP) for / and the points in An- 
In one variable, VDM{Ani, AnN) 7^ provided the points in An 
are distinct. Given a compact set K C C"', we say that K is determin- 
ing for [jVn if whenever h E [jVn satisfies h = on K, it follows that 
h = 0. For these sets we can find points {Ani, AnN} for each n with 
VDM{Ani, AnN) 7^ 0; we call these points unisolvent of degree n. De- 
spite the lack of a Hermite-type remainder formula, we can describe one 
condition on an array {Anj}j=i^...^N; n=i,2,... lying in a compact set K d 'C'^ 
satisfying a multivariate version of "regular with C \ connected" which 
implies for any / which is holomorphic in a neighborhood of K, we have 
Lnf ^ f onK. We call 

N 

An := supA„(2;) := sup |/„j(2;)| 

the ra-th Lebesgue constant for K,An (the function z — ?■ A„(z) is the n—th 
Lebesgue function). It is the norm of the hnear operator £„ : C{K) Vn C 
C{K) where Cn{f) = Lnf from (11 ■4p and we equip C{K) with the uniform 
norm. The next result follows from a multivariate version of Theorem [T] 
together with the Lebesgue inequality which says that for every continuous 
function f on K we have 

(1.5) \\f-Lnf\\K<{l + An)dn{f,K) 

with dn{f, K) as in fll.2p using polynomials in C^. 

Proposition 3. Let K G be polynomially convex and L— regular and let 
AnC K satisfy VDM{Ani, A„Ar) ^0 for each n = 1,2, ... //lim SUPyj_^QQ An 
1, for each f G C{K), 

hmsup 11/ - L„/||^/" = limsuprf„(/,ir)i/". 

n— >oo n— >oo 

The notions of polynomial convexity and L— regularity will be defined in 
the next section. 

This property lim SUp„_^QQ An — 1 is one of several we consider in the 
definition below. 

Definition 4. Let K be compact. Consider the following four properties 
which an array {Anj}j=i,...,Af; n=i,2,... C K may or may not possess: 

(1) lim^^oo A^" = 1; 

(2) limn^^\VDM{Ani,...,AnN)\'^''^ = 6{K); 

(3) lim^^oo jf EJLi ^a^j = fJ'K weak-*; 

(4) Lnf ^ f on K for each / holomorphic on a neighborhood of K. 
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Here 5a denotes the unit point mass at A. The probabihty measure 
liK is the (pluri-) potential theoretic equihbrium measure of K] i.e., for 
K d nonpolar, = ^'^^k K G C'^ nonpluripolar with d > 1, 

fj,K = j^^^idd'^VRY, the complex Monge-Ampere measure of V^{z) := 
limsup^_^2 Vft:(C) where 

Vk{z) = sup{^^log|p(z)| :pe\JVn, WpWk < 1} 

= sup{u(z) : u G L(C^) and M < on K}. 

Here L(C'^) is the set of all plurisubharmonic functions on C"' of logarithmic 
growth; i.e., u G L{C^) if u is plurisubharmonic in and u{z) = log \z\ + 
0(1) as l^l — )■ oo. In [8] it was shown that for i^" C C regular; i.e., Vk = 
= qk, and C\K connected, we have the implications 

(1.6) (1) =^ (2) =^ (3) =^ (4), 

while none of the reverse implications are necessarily true (although for 
arrays lying on the boundary of K, (3) and (4) are equivalent; for a more 
precise discussion, see [H]). Proposition |3] being true in for any d shows 
that the implication (1) =^ (4) remains true in C''; and, as was shown in 
[8], (1) =^ (2) as well. 

We continue in the next section with the necessary definitions and an 
elaboration on the relationship between conditions (1) and (2). Recent deep 
results of R. Herman and S. Boucksom ([3j and with Nystrom |3]) yield 
that (2) =^ (3); we discuss consequences of this result on recovering the 
measure fix in section 3. In section 4 we describe methods of recovering the 
extremal function Vk- We discuss the important Hernstein-Markov property 
in section 5. A brief introduction to weighted pluripotential theory in 
is provided in section 6, and a connection with unweighted pluripotential 
theory in C''^^ as in [7] is given. Section 7 provides explicit and semi-explicit 
constructions of arrays in certain compact sets satisfying conditions related 
to (l)-(4). We give a reprise of the analysis of so-called Bos arrays on the real 
unit disk i?2 C C in section 8. In section 9, we discuss computational 
approaches to constructing arrays in a compact set K satisfying some of 
the properties (l)-(4). A brief discussion of Kergin interpolation forms the 
content of section 10, and we conclude this work, as was done in [8], with a 
list of ten open problems. 

We would like to thank the organizers of the Conference on Several 
Complex Variables on the occasion of Professor Jozef Siciak's 80th birthday 
for their hospitality and we dedicate this work to Professor Siciak for his 
contributions and inspiration to the pluripotential theory community. 

2. SUBEXPONENTIAL LEBESGUE CONSTANTS AND ASYMPTOTIC FEKETE 

ARRAYS. 

We work in using the same notation as in section 1. For a compact 
subset K CC^ let 



Vn = Vn{K) := max \VDM{Ci, ...,Cn)\. 
Ci)---)CjvGX 
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Then 

(2.1) 6{K) = 6\K) = hm K 



dnN 

n 



is the transfinite diameter of K. Points zi,...,zn € K satisfying Vn = 
\VDM{zi, ...,Zf^)\ are called n—th order Fekete points for K. The tempo- 
rary superscript "1" refers to a weight w = 1 (see section 6). Zaharjuta [42] 
showed that the limit exists. Clearly if a compact set K is contained in an 
algebraic subvariety of C"' then 5{K) = 0. It turns out that for K G 'C'^ 
compact, 5{K) = if and only if K is pluripolar [,29J. 

If the compact set C C"' is L-regular, meaning that Vk = V^, and for 
i? > 1 we define 

(2.2) Dn:= {z:Vk{z)< log R}- 

then we have the Bernstein- Walsh inequality 

\p{z)\<\\p\\kR'"'^^, zEDr 

for every polynomial p in C^. A compact set C C"' is polynomially convex 
if K coincides with its polynomial hull 

K := {z G C''' : \p{z)\ < \ \p\\k, P polynomial}. 
Then Theorem [1] of Walsh goes over exactly to several complex variables: 

Theorem 5 ([36]). Let K be an L-regular, polynomially convex compact set 
in C^. Let R > 1, and let Dr he defined by 1^2. Let f be continuous on 
K. Then 

limsuprf„(/,is:)i/" < l/R 

n— >oo 

if and only if f is the restriction to K of a function holomorphic in Dr. 

The "only if" direction is the same for any d and uses the Bernstein- 
Walsh inequality. Theorem [5] immediately yields Proposition [3], showing 
that for L-regular, polynomially convex compact sets in C^, condition (1) 
implies (4). Unless otherwise noted, when discussing conditions (l)-(4) of 
Definition m we will always assume K is L-regular and polynomially convex. 

It is easy to see that (1) implies (2) but the converse is not true. On pp. 
462-463 in [8], it was observed that for an array {v4„j} C K with 

\VDM{Anl,...,AnN)\=CnVn{K) 

where 

< c„ < 1, limsupcy™ < 1, and lim c^'" = 1 

(e.g., c„ = for < f < 1), property (2) holds but (1) does not. More 
precisely, we have the following. 

Proposition 6. Let {Anj}j=i,...,Ar; n=i,2,... d K be an array of points. Sup- 
pose that 



,. t Vn{K) Y''' = \ 



Then (1) holds. 
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Proof. The result follows trivially from the observation that if 

\VDM{A^u...AnN)\ - 

then An < N ■ an- This observation is a consequence of the fact that each 
FLIP can be written as 

VDM{AnU...,Z,...,AnN) 
VDM{AnU...,AnN) 

SO that 

\VDM{Anl,...,Z,...,AnN)\ 



\Lj{z)\<an ^^^^^ 
Since \ VDM{Ani, z, AnN)\ < Vn{K) for each z E K,we have | |/rij| |x < 



3. Arrays yielding fix- 

In one variable, ^\og\VDM{zi, Zn)\ is a discrete "approximation" 
to the logarithmic energy of the measure /i„ = ^ X] j'=i ■ This is the idea 
behind the classical proof that (2) =^ (3). In several complex variables, 
the complex Monge-Ampere operator is non-linear and, until recently, no 
reasonable notion of the energy of a measure existed. We state without 
proof the remarkable result of Berman, Boucksom and Nystrom [1] that, 
nevertheless, (2) =^ (3) for general nonpluripolar compact sets K in the 
multivariate setting. 

Theorem 7 Let K G be compact and nonpluripolar. For each n, 

take points x^['\x^\ ■ ■ ■ , x^-* G K for which 

(3.1) lim \VDM{x^;'\--- = 6{K) 

rt— >oo 

(asymptotic Fekete points) and let /i„, := jfYlf=i^ (")■ Then 

fin — fiK weak — *. 

This gives a positive answer to question 6 posed in [8]. In Proposition 
3.7 of [8] it was shown that for a Leja sequence {xi,X2, ...} C K, 

lim \VDM{xu-,xn)\^ = S{K). 

n—>oo 

Thus the asymptotic Fekete property (13. ip holds for this sequence of points; 
so from Theorem [7] it follows that the discrete measures 

1 ^ 

i=l 

converge weak-* to fix- A Leja sequence is defined inductively as follows. 
Take the standard monomial basis {ei,e2,...} for [jVn ordered so that 
degCj <degej if ^ < j. Given m points zi, Zm in C^, we write 

VDM{Zi,...,Zn,) = det[ei{Zj)]ij=i^,„^rn- 
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Starting with any point Xi G K, having chosen Xi,...,Xm G K we choose 
Xm+i & K SO that 

\VDM{xi, ...,Xm,Xm+i)\ = mayi\VDM{xi, ...,Xm,x)\. 

We remark that despite possessing the desirable property that fin Hk 
weak-*, it is unknown if (1) always holds for a Leja sequence, even in the 
univariate case {d = 1). This is the first question in section 5 of [8]. We end 
this subsection with the statement of a result of R. Taylor and V. Totik 
that gives a partial answer in the d = 1 setting. 

Theorem 8 ([39j). Let K G C be compact and assume that the outer 
boundary of K can be written as a finite union of arcs. Then any Leja 
sequence for K satisfies property (1). 

In particular, Leja sequences on an interval satisfy property (1). We 
return to this topic in section 7. 

4. Recovering the function Vk- 

In section 2.4 of [8], an elementary argument showed that for arrays 
satisfying property (1) of Definition HI the Lebesgue functions A„(z) can be 
used to recover Vk in the sense that 

(4.1) lim logA„(z) = Vk{z), z e C. 

This was proved in the univariate case but the same proof works in all 
dimensions. Property (2) is not sufficient for (14. ip to hold. In this section, 
we investigate special families of polynomials which can be used to recover 
the extremal function Vk- 

Zaharjuta's proof of the existence of the limit in (12. ip introduced the 
useful notion of directional Chebyshev constants. Let ei{z), ...,ej{z), ... be a 
listing of the monomials {ei{z) = z"^*^ = • • ■ z^''-} in indexed using 
a lexicographic ordering on the multiindices a = a{i) = (ai, G N'^, 

but with degCj = \<y{i)\ nondecreasing. Define the class of polynomials 

P, = P(a(0):={e,(^) + ^c,e,(;.)}; 

j<i 

and the Chebyshev constants 

T{a) ■.= mi{\\p\\K:peP.}. 

We write ta,K '■= ioi(i),K for a Chebyshev polynomial; i.e., t^^K ^ Pi and 
= T{q.). Let S = denote the standard simplex in W^] i.e., 

d 

T. = {9 = (^1, ...,9d) -.Y, 9, = 1, 9^ > 0, J = 1, 4, 

i=i 

and let 

S° := {0 G E : 9^ > 0, j = l,...,d}. 
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For all 6* e S°, the limit 

T{K,e) := lim T(a)^/I"l 

a/\a\—>-0 

exists and is called the directional Chebyshev constant for K in the direction 
9. Zaharjuta showed that 

(4.2) 5{K) = exp [ \— [ log r (iT, e)d9] . 

'-meas(Lj Jj^o 

In [B], the following theorem was proved. 

Theorem 9 ([6J). Let K G be compact, L— regular, and polynomially 
convex. Let {pj} be a sequence of polynomials such that for all 6 G S*^, there 
exists a subsequence Yg C Z+ with pj e -P(«j), j G Yg and 



Then 



lim|b,|ir^^^=r(ir, 



[linisup-^log^M]* = K^(^), ,^K. 



The family {pj} is said to be 6 — aT - "theta-asymptotically Chebyshev" 
- if the property in Theorem |9] holds. Bloom proved, in particular, that Leja 
polynomials associated to a Leja sequence have this 6 — aT property. Using 
Theorem ini he proved an interesting corollary related to our condition (1). 
To this end, we begin with a triangular array {-Bsj}j=i,...,s; s=i,2,... with the 
property that VDM[Bsi, Bss) ^ for each s. Define, for each multiindex 
a = a{s), the polynomial 

VDM{Bsu...,Bss,z) 



VDM{Bsu-.,B,s) ■ 

Note that Ga(s){z) = z"'^^') + ^j^gCjz"^^^ G Ps- Moreover, it is straightfor- 
ward to see that 

Ga{s){z) = Z"^""^ — La{s){z°'^''^) = ta{s),K " L a{s){t a{s) ,k) 

where La(s){f) is the LIP for / and the points -B^i, ...,Bss; i.e., 

s 

L^^,){f){z) = Y,f{Bs,)lsj{z) 



and 



Letting 



VDM{Bsi, -Bs(j-i), z, .--Bs 
''^'^ " VDM{B,,,...,B,.,) 



Aa{s) ■= SUp^ 

we have the following. 
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Theorem 10 ([6J). If K is L— regular and lim\a{s)\^oo -^^^^^^^'^'^^ = 1; then 



Vk{z) = limsup r^-^ log 77— r— 

|a(s)|^oo \Oi[S)\ \\<^a{s)\\K 



for zeC^XK. 



The hypothesis hm|Q,(s)|_>oo ^1(1')''^^' ~ -'- ™pli6s that the family of polyno- 
mials {Ga{s)} is 9 — aT. In particular, Fekete polynomials for each s = 1, 2, ... 
defined from an array that maximizes \ VDM{(i, ...,Cs) \ over {(i, ...,(s) € 
are shown to have this property. Note that in this case, for s = N = dimP„, 
the s points Bgi, Bgg coincide with the n— th order (degree) Fekete points 
AnN- A weighted version of Theorem [9] was proved as Theorem 3.5 
of [13]. We will utilize this in section 6. 



5. Bernstein-Markov property. 

For a compact set K G and a measure u on K, we say that the pair 
{K, u) satisfies a Bernstein-Markov property if there exist constants {M„} 
with lim sup^^oQ M^" = 1 and all polynomials Qn € Vn satisfy 

n 1 1 L'^iv) • 

In [8] it was shown (Theorem 3.3) how one could recover the transfinite 
diameter 5{K) from asymptotics of Gram determinants associated to a 
Bernstein-Markov pair {K,v). More recently, strong Bergman asymptotics 
were proved in [3] in this setting: if [K, v) satisfies a Bernstein-Markov 
property, then 

(5.1) —B^dv — weak-* 

where 

TV 

i=i 

is the n— th Bergman function for v and {gi, 5'2, ■ ' ' i Q'Af} is an orthonor- 
mal basis for Vn with respect to L'^iy)- Thus it is natural to ask which 
compact sets K admit measures v satisfying a Bernstein-Markov property. 
The following result was proved in [15]; since the proof is short, we include 
it. 

Proposition 11 ([15]). Let K G be an arbitrary compact set. Then there 
exists a probability measure v such that {K, u) satisfies a Bernstein- Markov 
property. 

Proof. To construct we first observe that if is a finite set, any measure 
V which puts positive mass at each point of K will work. If K has infinitely 
many points, for each k = 1,2,... let =dimVk{K), the holomorphic 
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polynomials on restricted to K. Then linifc^^oo f^k = oo and rrik < {'^^^) 
0{{d)''). For each k, let 

ink 

E5 (fc) 
2 



where {2;]''''}j=i,...,mfe is a set of Fekete points of order k for K; i.e., if 
{ei, e^fe} is any basis for Vk{K), 

(5.2) |det[ei(z-^^)]ij=i,„,,mJ = max |det[ei(gj)]i,i=i,...,r; 



<Jl,...,IJmi.6^ 



Define 

where c > is chosen so that u is a probability measure. If p G Vk{K), we 
have 

i=i 

where G Pfc(ir) with /f^l^''^) = ^jfc- We have ||/f^||i^ = 1 from flOjl 
and hence 

lblk<Elp(#^)l- 

On the other hand, 



17/1/'^ 

krukilogky ^ 



Thus we have 



A;mfc(log/i;)2 

blk < m\L\du)- 



□ 



This gives a positive answer to the first part of question 2 in [8J. The second 
part has a negative answer, as the simple example of 

K = {z eC:\z\<l}VJ {2} 

shows. 

For certain measures v with compact and non-polar support on the real 
line M C C, pointwise asymptotics of the Bergman functions {-B^} are known 
(cf., [3Q]). In the higher dimensional setting, very little is known. Bos, et al 
[T8] consider one natural analogue of the interval, namely, the real unit ball 

d 

Bd := {(^1, Zd) e C : Im;^,- = 0, j = 1, d- J](Re2,)' < 1} 
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in M*^ C C''. Writing Xj := Rezj and x = (xi, ...,Xd), it is known that 

2 

d^Bd = ^ {x)dx := dx 

where dx = dxi A ■ ■ ■ A dxd is d— dimensional Lebesgue measure and Ud is 
the surface area of the unit sphere S'^ C M'^"'"^. Lemma 1 in [T8] shows that 
if du^x) = u{x)dx where u{x) = u{—x) is a positive centrally symmetric 
weight satisfying a certain Lipschitz property, then 

(5.3) hm i-i^r^-^ 



and the convergence is uniform on compact subsets of the M*^— interior of 
Bd- The first step is the special case where u{x) = (note in this case 

the right-hand-side of (15. 3p is the constant function 1; it is also shown that 
the limit is the constant function 2 on the sphere |x| = 1). 

Having a Bernstein-Markov measure allows one to replace Chebyshev 
polynomials by orthogonal polynomials in certain asymptotic computations. 
Using this idea, the exact calculation of the transfinite diameter of the real 
ball Bd := {x E M."^ : |x| < 1} and the real unit simplex, S = Sd '■= {x E 
^+ ■ Yl'j=i^j — ^} ^ '^'^ recently achieved in |9]. 

Proposition 12 ([9J). The transfinite diameter of the unit ball Bd is: 
(1) for d even, 

1 / 12rf + lAl 1 1, I 



wr, N 1 / 12rf+l v^l 1 1, 1 
«™ = 5exp -J— ^-H-- + -log(2) + -5; . 

\ i=i i=i 

(2) for d odd, 

1 / 12rf+l>Al 1 d-1 iJ^i-lV 



, 1 / 12rf+ 1 1 1 d - 1 , , 1 
5{Bd) = -exp ---- - + - + ^log(2) - - 



2--^l 4 d ^j^2^ 2d ^^^^^^ Ad^ J 

The transfinite diameter of the simplex Sd is 

8{Sd) = {S{Bd)f. 

6. Weighted vs. unweighted. 

In the weighted theory, one considers closed sets which, for certain weights, 
may be unbounded. To be precise, let K C he closed and let w be an 
admissible weight function on i^: w is a nonnegative, use function with 
{z E K : w{z) > 0} nonpluripolar; if K is unbounded, we require that w 
satisfies the growth property 

(6.1) |z|w(z) — )• as |2;| — )■ oo, z E K. 

Let Q := ~ log w and define the weighted extremal function or weighted 
pluricomplex Green function V^q{z) := \\msvL]i(-^^VK,Q{C) where 

Vk,q{z) := snp{u{z) : u E LiC^), u < Q on K}. 
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In the unbounded case, property (16.11) is equivalent to 

Q{z) — log \z\ — 7- +00 as \z\ — 7- oo through points in K. 

Due to this growth assumption for Q, Vk,q is well-defined and equals VKnBR,Q 
for R > sufficiently large where Br = {z : \z\ < R} (Definition 2.1 and 
Lemma 2.2 of Appendix B in [33]). It is known that the support 

Su, := supp(/ix,Q) 
of the weighted extremal measure 

1 



is compact; 



moreover, 



S,^CS*^:={zEK: VIq{z) > Q(z)}; 



Vk,q = Q q-e- on 5*^ 
(i.e., V^Q = Q on Sw\F where F is pluripolar); and if m G L(C'^) satisfies 
u < Q q.e. on Sw then u < V^ q on C^. 

The unweighted case is when K is compact and w = 1 {Q = 0); we then 
write Vk '■= Vka to be consistent with the previous notation. 

Even in one variable {d = 1) the weighted theory introduces new phe- 
nomena from the unweighted case. As an elementary example, puts no 
mass on the interior of K (in one variable, the support of fix is the outer 
boundary of K); but this is not necessarily true in the weighted setting. 
As a simple but illustrative example, taking K to be the closed unit ball 
{z : \z\ < 1} and Q{z) = it is easy to see that Vk,q = Q on the ball 
{z : 1^1 < 1/^2} and VK,Qiz) = log \z\ + 1/2 - log(l/\/I) outside this baU. 
One can check that if is the ball {z : \z\ < R} for 1 < R < oo and 
Q{z) = one has Vr-^.q = Vk,q. 

Now let i^' C be compact and let w be an admissible weight function 
on K. Generalizing (11. 3p . given d, ...,(n £ K, let 

W{Cu Cn) := VDMiC, CN)w{C^r ■ ■ ■ HCnT 

ei(Ci) ei(C2) ••• ei(Civ) 



det 



wiCiT-'-wiCNY 



. e7v(Ci) eAr(C2) • • • eiv(CAr) . 

be a (generalized) weighted Vandermonde determinant of order n. Let 

W^{K) := max |W^(Ci,...,Civ)| 
Ci,...,Cjve-ff 

and define an n—th order weighted Fekete set for K and w to be a set of 
points Ci, Cat G K with the property that 

\W{Ci....XN)\ = Wn{K). 
We also write (5"''"(i^) := Wn{K)^ and define 



1 dnN 
■ ' n\^^ J 

n— ^-oo n— >oo 



5"'(ir) := lim 6'"^''{K) = lim Wn{K) 
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A proof of the existence of the hmit may be found in [3] or [T3]; in the latter 
work one defines the circled set 

F = F{K, w) := {{t, z) = {t, tX) e C'^+^ -.XeK, \t\= w{X)} 
and shows that, indeed, for the closure F, 

6'"{K) = 6(F)^. 
In |1] the authors proved a weighted version of (2) =^ (3): 
Theorem 13 (P]). Let K G C'^ be compact with admissible weight w. For 
each n, take points x^i\x^2 \ ' ' ■> ^n^ ^ ^ /^'^ which 

lim \W{x^;'\--- ,x^;;^)\^ =6^{K) 

n— >oo 

(asymptotically weighted Fekete points) and let /i„ := J2f=i ^^f")- Then 

fJ-a IJ'K,Q weak - *. 

The main results in [3], which are stated and proved in a much more 
general setting than weighted pluripotential theory in and lead to the 
results in require the weighted theory. A self-contained exposition of the 
weighted pluripotential theoretic setting can be found in |28] . 

As a final remark on weighted pluripotential theory, we provide a solution 
of Problem 3.4 in Appendix B of |33j in the locally regular, w continuous 
case, correcting the end of section 8 in [13]. The problem is to show that 
the weighted Fekete polynomials can be used to recover the weighted pluri- 
complex Green function Vk,q in the sense of Theorem 3.5 of [13]; see f l6.2p 
below. Let i^' C be locally regular and let w be a continuous admissible 
weight on K. We define weighted Chebyshev constants 

T^{K) := M{\\w\''^%\\k : p G P^}^ 
and we let t^ denote a weighted Tchebyshev polynomial; i.e., is of 
the form w'^'-'^p with p e Pi and \ K^k\\k = rf(/s:) 

Let nig = ('^^*). Note that \a{i)\ = s for m^.i + 1 < i < m^. Given 
Ci,...,0 e K, let 

WiC, C.) := VDM{(,, O)^(Ci)'"^'^' ■ ■ ■ ^(0)'"^'^'. 

Generalizing the notion of an n—th order weighted Fekete set for degree n, 
for each i = 1,2,..., an i—th weighted Fekete set for K and w will be a set 
of i points (i, ^ K with the property that 

Wr.= \W{Cu...,0)\= sup 

Fix i with \a{i)\ = s. We will define weighted Fekete polynomials pi for each 
positive integer i with \a{i)\ = \a{i — 1)| = s; i.e., rng-i + 1 < « — 1 < mg. 
Choose an (z — 1)— st weighted Fekete set (i, O-i and form the weighted 
polynomial 

w{zr^%,{z) = w{zr^^He,{z) + J2c,e,{z)} : ^(Ci' O-i, ^) 



j<i *■ 
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Thus 

Next we choose an i—th weighted Fekete set In the expansion of 

the determinant of this weighted Vandermonde 

W{C,, Q) := VDMiC^, C.)^(Ci)'"^^^' ■ ■ ■ ^(0)'"^^^ 

we replace the last row w(Cfc)'"^*-"ei(Cfe), k = l,...,i by t^(j)^;^(Cfc), k = 1, i. 
Expanding the weighted Vandermonde determinant by the last row, 



W,<J2Kir),Km-\W{Cu...,Q-i\ 
1 

<irf(i^^)l"»liyi_i. 



k=l 



Thus 



The sequence of weighted polynomials {w'^^^^lpj : a{i) ^ (s,0, ...,0), s = 
1, 2, ...} satisfies the hypothesis of Theorem 3.5 of [13j - note that (1, 0, 0) ^ 
E° - so that we obtain its conclusion: 

(6-2) [limsup^log '^l^f ]* = Vk,q{z), z^K. 

L i^oo |a(^)| \W'^^'>\pi\\K 

7. Explicit good interpolation points 

An important question in numerical analysis and computational math- 
ematics is to provide explicit or computable arrays {^4^} — {A^j, j — 
1, . . . , A^} C K satisfying (4) in Definition |H By this we mean points whose 
values can be entered on a computer with arbitrary precision to evaluate 
LIP's. Typical explicit points in the univariate setting are Chebyshev points 
for which 

{A„} = (cos f J^^vr^ , k = 0,...n], K=[-l,l], 



or the roots of unity for which 

{An} = {exp{2ikn/{n + l)), k = 0,...,n}, K = {\z\ < 1}. 

Semi-explicit points, such as Fejer points (the image of a complete set of 
roots of unity under an exterior conformal mapping) or points related to 
orthogonal polynomials still deserve interest since, in particular cases and 
for n not too large, they can be computed with high precision. The Fekete 
points for [—1,1] are extreme points of Legendre polynomials plus {—1,1}. 
Recently Xu [H] constructed semi-explicit points on a region in bounded 
by two lines and a parabola using zeros of Jacobi polynomials. Typical non- 
explicit good points are multivariate Fekete points, which cannot currently 
be efficiently computed. Even in the one variable case, only a few examples 
of explicit good interpolation points are known. If K is not an interval or a 
disk one generally uses algorithms which provide numerical approximations 
for the points (cf., section 9). 



POLYNOMIAL INTERPOLATION AND APPROXIMATION IN C* 



15 



To clarify what we mean by good points for polynomial interpolation we 
utilize the following hierarchy which slightly refines some conditions listed 
in Definition HI 

Definition 14. Let {An} be an array of interpolation points in K G C^. 
We denote by A„ the Lebesgue constant for An, K and L„/ the LIP of /, An- 
Here are four properties that {An} may or may not possess: 

(HI) Lnf ^ f on K for each / holomorphic on a neighborhood of K, see 
Definition m (4). 

(H2) The Lebesgue constants A„ grow sub-exponentially: lim„_>oo A.^" = 

1, see Definition m (1). 
(H3) The Lebesgue constants A„ grow polynomially: there exists s G N 

such that An = Oiji") as n — )■ oo. 
(H4) The Lebesgue constants A„ grow sub-polynomially: A„ = o(n) as 

n — )• oo. 

Clearly (H4) =^ (H3) =^ (H2) =^ (HI). None of the reverse 
implications is true. That there are sequences of points satisfying (H3) but 
not (H4) and (H2) but not (H3) will follow from the results below. In the 
multivariate case, since there is currently no analogue to Theorem [2], we 
are left with the last three conditions and are obliged to study multivariate 
Lebesgue constants. In the univariate case, we may establish (HI) without 
having recourse to Lebesgue constants, e.g., in obtaining a discretization of 
the equilibrium measure (using (3) =^ (4) in flL6p ). An array of multi- 
variate Fekete points is a fundamental example satisfying (H3). 

For theoretical approximation of holomorphic functions on a neighbor- 
hood of K, conditions (H4) and (H3) do not provide better results than 
(H2). However, computations with arrays of points having smaller Lebesgue 
constants benefit from a higher stability and, from the point of view of ap- 
proximation theory, one can derive convergence results for larger classes of 
functions. For instance, from Lebesgue's inequality (II. Sp and Jackson's the- 
orem, see e.g. |31j §1.1.2], if if is a product of intervals in C then 
(H3) implies that L„/ ^ f on K for each / which is (s -|- l)-times con- 
tinuously differentiable on a neighborhood of K while (H4) only requires / 
to be continuously differentiable. In general, conditions (H3) and (H4) im- 
ply convergence results for classes of functions / for which dn{f,K) (recall 
( II. 2p ) is known to decrease polynomially in n. Such estimates are known for 
several natural spaces of functions holomorphic on the interior of K with 
some regularity up to the boundary (cf., [34J and the references therein). 

The array of Padua points {Pn} is an explicit example of multivariate 
interpolation points satisfying condition (H4). We follow the presentation 
given in [19]. Another point of view can be found in [20]. The points of 
Pn lie in K = [—1, 1]^ and are located on the classical Lissajous curve 
t G [0,7r] I— 7- 7„(t) := (cos(nt), cos((n + l)t)). The Padua points are the 
double points of this curve together with its points on the (real) boundary 
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of the square. A simple formula is the following: 

(7.1) p^^= LJ^^n + ^Ti] : i + j<n, {i,j)eN' 
I \n + 1 n J 

It is readily seen that the above ("2 ^) points are pairwise distinct but the 
fact that P„ is a unisolvent set of degree n is not immediate. To see this, 
it suffices to exhibit a FLIP for any point a G P^- There is a remarkable 
formula expressing such a FLIP with the help of the reproducing kernel Kn 
for the inner product based on the tensor product of two arcsine measures, 



{p,q) eVn^ {p,q) ■= ^ P{.x,y)q{x,y) dxdy 



— ^ / ]9(cos(t), cos(s))g(cos(t), cos(s))(it(is. 



./[0,7rl2 



The reproducing kernel Kn is defined on [—1,1]^ via the relation 

It can be shown that the FLIP corresponding to a = {xa,ya) £ Pn is given 
by 

laix, y) = Wa {Knia; (x, y)) - Tn{y)Tn{ya)} 

where Wa is constant and T„ is the ordinary Chebyshev polynomial of degree 
n. The proof uses certain quadrature formulas for the the tensor product of 
two arcsine measures using the points of P„. Next, it can be shown that the 
kernel Kn, hence the FLIP's for P„, are expressible as a linear combination 
of quotients of classical trigonometric polynomials. A careful analysis leads 
to the following result. 

Theorem 15 ([19]). The Lebesgue constants A„ for the Padua points P„ 
and K = [—1,1]^ satisfy 

An = 0(log^ n), n ^ oo. 

In particular, Lnf ^ f on K for each f which is continuously differentiable 
on a neighborhood of K. 



Lagrange interpolants at Padua points can be easily computed, cf., 
. Unfortunately, it is not clear what might be the analogues of Padua 
points in higher dimensions (n > 2). 

We present another construction of good points based on a different 
idea. These points will only satisfy the weaker condition (113) but there is a 
simple and efficient way of going from dimension 1 to dimension k for every 
k. The starting point is a classical algebraic formula giving multivariate 
interpolation points starting with univariate points. Given d sets of (n + 1)- 
tuples 

As = (Oos, • • • , Clns) , S = 1, . . . ,d 

consisting of distinct points in C, we intertwine these tuples: 

Ai © Aa © ■ ■ ■ © Ad := {{ai^i, ai^2, ■ ■ ■ , ai^d) , «i + «2 H \- id < d} . 
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Changing the ordering of the points would provide a different set of points 
in C^. The set of points we obtain is unisolvent of degree n in C'*. Such 
points were studied by Siciak and used in [SJ Theorem 4.8] to prove 
the existence of points satisfying properties (2), (3) and (4) in Definition 
m We want to construct interpolation points by intertwining well-chosen - 
and well-ordered - univariate interpolation points which satisfy (H3). The 
obvious strategy is to try to relate the Lebesgue constants of the multivariate 
interpolations points to the Lebesgue constants of the univariate points. 
Such a relation is given in the following theorem. For simplicity, we state 
only the case d = 2. 

Theorem 16 ([23]). Let K be a compact set in containing Ai © A2. 
We let Ki {K2) denote the projection of K on the zi (22) axis. Then 

where Ai^'^ = (oq, . . . , Oj) and A„ (^Ai^^^\Ki^ denotes its Lebesgue constant 
with respect to the compact set Ki {likewise for A2'"''). 

In order to estimate the Lebesgue constant A„(Ai © A2\K), bounds on 
A„(Ai|_ft'i) and on A„(A2|/^2) do not suffice. We must have bounds on the 
Lebesgue constants of every subset Ai^*' and A2'''^ for i + j < n. The only 
practical way of using the theorem is to start with univariate points given 
by a sequence of interpolation points A" = (aos, . . . , ctns), so that for every 
i < n, Ag = Al. Then the search for good multivariate interpolation points 
via the intertwining process is reduced to the problem of finding univariate 
interpolation points given by a sequence and satisfying (H3). Surprisingly, 
such sequences did not seem to be known until recently. 

We now discuss the construction of such univariate sequences. All exam- 
ples currently available are constructed with the help of Leja sequences for 
the closed unit disk. Recall that a Leja sequence for a compact set K G C 
is a sequence {a„} in K such that 

d d 

max If \z — aA = |T Ifld+i — aA, d>0. 

i=0 1=0 

If we are to produce explicit points we must restrict to K = D = {\z\ < 1}. 
In this case, the structure of Leja sequences is given by the following result. 
We always assume that the first term equals 1. 

Theorem 17 (^Sj). Leja sequences for the unit disk D satisfy: 

(1) A 2^-Leja section is formed by the 2"'-th roots of unity. 

(2) //£'2n+i is a 2"'^^-Leja section then there exist a 2^-th root p of —1 
and a 2^ -Leja section U2" such that £2^+1 = {E2« , pf/2"). 

Using this result, the following estimates were recently established. 

Theorem 18 ([25]). Let {cj} be a Leja sequence for D . As n ^ 00, A„ = 
O(nlogn) where A„ is the Lebesgue constant for {cq, . . . , e„_i}. 
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A similar estimate holds for the image of Leja sequences under external 
conformal mappings C \ D — )• C \ K for K sufficiently regular, e.g., K 
bounded by a Jordan curve It is also shown in |25] that a Leja 
sequence for D cannot satisfy (H4), thus showing that, in general, (H3) 
does not imply (H4). 

For practical applications, real points are more useful. A simple idea to 
construct such points is to project a Leja sequence for D onto the real axis. 
Since a Leja sequence for D is symmetric with respect the real axis, complex 
conjugate points provide the same real point. Eliminating this redundancy 
we obtain a so-called Re-Leja sequence [26]. One can specify the n-th entry 
of a Re-Leja sequence in terms of the real part of a certain entry of the Leja 
sequence used in its construction. 

Theorem 19 ([26j). Let X = {xj} be a Re-Leja sequence. The Lebesgue 
constants A„ for the points xq, . . . , Xn-i satisfy 

(7.2) A„ = 0(n^log n), n — i- oo. 

By combining the above two results and the general version of Theorem 
[16] we obtain the following result. 

Theorem 20 ([25l[26]). Intertwining k Leja sequences for D andd — k Re- 
Leja sequences for [—1, 1] yields unisolvent sets on the Cartesian product of 
these sets in whose Lebesgue constants A„ satisfy (H3). 

Moreover, the degree of the polynomial growth of A„ can be estimated. 
One can give an explicit expression for the n-th element of certain simple 
Leja (or Re-Leja) sequences that depend only on the binary expansion of 
the index n. Details can be found in [251 12S] • In view of Theorems (TB] and 
[HI the intertwining of Leja sequences for many compact sets provides fur- 
ther examples of multivariate interpolation points satisfying property (112). 
Goncharov [27] has constructed a sequence in [—1, 1] satisfying (H2) but 
not (H3) by arranging the classical Chebyshev points in a certain manner. 

8. The quest for good points in the real disk 

We describe a natural strategy for finding good points in the real disk 
B = B2 = {{x,y) G : x'^-\-y'^ < 1}. Using Theorem 3.3, one can calculate 
the transfinite diameter of B to find 

m - ^ 



'2e 

and its equilibrium measure is well-known, 

r 

dfiB = - — ^ drdO 



r 



in polar coordinates in M , i.e., d^iB is absolutely continuous with respect 
to Lebesgue measure on B with density 



It was shown in [8] that (3) ^ (1), (2), (4). The set B was used in con- 
struction of Bos arrays satisfying (3) but not (4) (and (3) but not (2)). The 
points at the n— th stage in a Bos array are formed by taking a union of 
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equally spaced points on concentric circles centered at the origin. Precisely, 
if n = 2s is even, one chooses s + 1 radii RgQ < ■ ■ ■ < Rss = 1 and 4j + 1 
equally spaced points on the circle of radius Rsj. The Vandermonde deter- 
minant \VDM{Ani, Aniy)\ dcpeuds only on the radii Rso,--- ,Rss and 
if the asymptotic distribution of the radii on [0, 1] is given by a function 
G : [0, 1] ^ [0, 1]; i.e., if G'(^) = R%, then 

lim \VDM{Ani,...,Ar,N)\'^'" = -^exp 



V2 ' AL{G) 
where 

L{G) = I x^\ogG{x)dx + 2 f f x\og[G{y) - G{x)]dydx. 

Jo Jo Jx 

Thus, if one could construct a Bos array with L{G) = —2/3, then this array 
would satisfy (2). 

Taking G{x) = (1 — cos7ra;)/2, the radii distribute asymptotically like 
the Chebyshev distribution on [0, 1] and this is a necessary condition (see 
[S]) that such an array satisfies (4). We state without proof an interesting 
calculation. 



Lemma 21. For G{x) = (1 — cos(7rs))/2 = sin^(7rx/2) 

4 2 

-7Tlog(2) + - 
O TT 

where ({x) is the classical zeta function. 



L{G) = -- log(2) + 3^C(3) ^ -0.6806085842 • • • 



In particular with this G, L{G) ^ —2/3 so such a Bos array does not satisfy 
(2). 

Taking G{x) = 1 — (x^ — 1)^, the arrays satisfy (3): we obtain fiB as the 
limiting measure. Elementary but nontrivial calculations yield 

Lemma 22. For G{x) = 1 - (x^ - 1)2, 
26 

L{G) = -— - 4 log(2) + 4^2 log(y2 + 1) ^ -0.675675691 ■ ■ • . 

Thus again, with this G, L{G) ^ —2/3 so such a Bos array does not 
satisfy (2). Indeed, if one could find a Bos array with L{G) = —2/3, then, 
by Theorem [71 the array satisfies (3) and hence, a posteriori, G{x) = 1 — 
(x^ — 1)2, a contradiction. Thus, unfortunately. Bos arrays on B never satisfy 
(2), giving a negative answer to question 9 of [8j. 

9. Algorithms 

It is clear that one can expect to have lists of explicit good interpolation 
points only for a very limited class of compact sets, even in the univariate 
case. If we have to produce good points for a more or less arbitrary compact 
set, one must produce them algorithmically. We now discuss some recent 
work in this direction. 

In a series of papers. Bos, Sommariva and Vianello (cf. [38j and with De 
Marchi [T7j) have introduced the notion of approximate Fekete points. For 
K G compact, a basis {Pi,P2, ■ ■ ■ ,Pn} ioi Vn, and a set of M > 
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points {oi, qm} of K, we consider the N x M matrix whose columns are 
of the form 



Selecting a subset of columns is then equivalent to selecting a subset of 
points. We choose the first point xi G {ai, ...,aM} to maximize ||^(aj)||2. 
Having chosen Xi,X2,--- ,Xk G {ai, ...,aM} the {k + l)st point Xk+i G 
{ai, om} is chosen so that the volume generated by the columns V{xk+i) 
and V{xi), V{x2), ■ ■ ■ , V{xk) is as large as possible. 

Suppose K is L— regular. If for each n = 1,2,... one chooses a set of 
M(n) > points AM{n) = {a["\ aj^|^)} of K so that IJn^M(n) forms 
a weakly admissible mesh for K (WAM), then the corresponding array of 
approximate Fekete points satisfies (2) and hence (3) (Theorem 1 of [TB]). 
The mesh ^M(n) is weakly admissible, according to [21], if #^Af(n) grows 
polynomially in n and 

M\K<CnMW,,^, for ^WpeVn 

where Cn grows polynomially in n. All L— regular compact sets K admit a 
weakly admissible mesh; cf.. Theorem 16 of [24] . We remark that a WAM 
is called admissible (AM) if one can take C„ = C, a constant independent 
of n. 

There is also an algorithmic notion of discrete Leja points; as with ap- 
proximate Fekete points, constructing discrete Leja points from a weakly 
admissible mesh Un^A^(n) gives an array satisfying (2) and hence (3). The 
interested reader is referred to [T^ for details of the algorithm. 



10. KERGIN INTERPOLATION 

Of the many polynomial interpolation alternatives to Lagrange interpola- 
tion, one of the most productive ones utilized for interpolating holomorphic 
functions in C^, ci > 1 is Kergin interpolation. 

In this section, we give a new presentation of Kergin interpolation which 
highlights its canonical character. We let 0{C^) denote the space of entire 
functions and £(C(C'^), P„) the space of continuous linear maps from 0{C^) 

toVn. 

Theorem 23. There exists a unique map K, 

(10.1) /C : A = (ao,...,aOG(CY+i^/C[A]G£(0(C"),P.), 
such that 

(Kl) for every f G O(C^), /C[A](/)(a,) = /(a,), j = 0, ...,n; 
(K2) for every f G 0{'C^), the map A — > /C[A](/) is continuous; 
(K3) K, is coordinate- free. 
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The map JC is defined by 

n „ / ^ \ 

(10.2) /C[A](/)(x) = ^ / D'f lj2uaA{x-ao,...,x-ak.,)dmk{t), 

k=0 \j=0 J 

where D^f is the k-th total derivative of f, = {t = (to, • • • , ^fc) ^ 
[0, 1]'^"'"-'^ : Yli=o^i ~ ^} ^^'^ '^^fc Lebesgue measure on Sk- 

That /C[A] is coordinate-free means that for every invertible hnear map 
m on 

/C[A](- o m) = /C[mA](-) o m, 

where mA := (m(ao), . . . , m{ad)). 

Condition (Kl) is relatively weak (e.g., if = a for z = 0, . . . , d there is 
only one condition). We shall see later that, together with {K2) and {K3), 
it implies much stronger properties. The operator /C[A] is called the Ker- 
gin interpolating operator with respect to A. In contrast with multivariate 
Lagrange interpolation, the number of points n + 1 is independent of the 
dimension c? of C''. 

Proof. We first prove that there exists at most one map 

n : A G (C'^)"+i ^ n[A] G £(C(C'^),P„) 

satisfying (Kl), {K2) and {K3). Suppose Hi and 112 are two such maps. 
We prove that for every A G f/ := {A G (C^)*^"^^ : ai ^ aj for / 7^ j} and 
every / G 0(C'^), ni[A](/) = n2[A](/). Since U is dense in (C'^)"+i and 
A — )■ nj[A](/) is continuous, this suffices to prove our claim. 

We can reduce the problem as follows. Since nj[A] is continuous on 
0{C^) and the space V spanned by ridge entire functions - functions of the 
form z ^ h{{X,z)), where h G C(C) and A G - is dense in 0(C°'), it 
suffices to prove that Hi [A] = 112 [A] on V. Further, since Ilj [A] is linear we 
simply need to prove ni[A](/) = n2[A](/) for / = h{{X, ■)), with h G C(C) 
and A G C^, A ^ 0. 

Fixing such an /, let if = {(A, ■) = 0} be the hyperplane orthogonal 
to A. For e > we define a linear map by me(A) = A and m^ln = eld 
where Id denotes the identity on H. Clearly is invertible. Moreover, since 
m^x — X G if, we have 

(/ o me)(x) = /i((A, m^x)) = h{{X, x)) = f{x). 
Since Ilj is coordinate-free, we deduce that 

naA](/) = naA](/om,) = nam,A](/) om„ e > 0. 
We have me{x) — )■ (A,a;)pjp- as e — 0, thus by (-^^2), 

A 



(10.3) n4A](/)(x) = n, 



(AA) 



lAI 



(/) (A,a:) 



lAI 



where (AA) ■ pip = ((Aaj) ■ pjp : z = 0, . . . , ra). Since n, takes values in Vm 
(110. 3p implies that there exists a univariate polynomial p of degree at most 
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n depending on /, A and A such that 

n,[A](/)(x)=p((A,x)). 

We speciahze to the case where the (A, a,) are distinct. Since the Oj them- 
selves are distinct, the set U of all such A is dense in C^. It remains to use 
assumption {Kl). We have 

/i((A,ai)) = /(oi) = ni[A](/)(ai) =p((A,ai)), i = 0,...,n. 

Hence p is a polynomial of degree at most n that interpolates h at these 
n + 1 points, i.e, p„ is the LIP of h at these points which we write as 

(10.4) naA](/) = L[(A, ao), ... (A, a„) ; /i]((A, ■)), \eU. 

In particular, ni[A](/) = n2[A](/). We now use the density of U and the 
continuity of / — )■ nj[A](/) to extend the identity to the case where X ^ U. 
This finishes the proof of the uniqueness. 

Identity fll0.4p shows that if a map 11 with the required properties exists 
then it should come as a natural multivariate generalization of one of the 
many available expressions of univariate Lagrange-Hermite interpolation. 
Formula (110. 2p is the natural multivariate version of the classical Hermite- 
Genocchi formula. The proof that this map satisfies the required properties 
is a simple calculation; cf., ^30j- CD 

It is not difficult to show that the map /C[A] interpolates in the Hermite 
sense; i.e., if a point a appears k times in A then }C[A]{f){a) = f{a), 
j = 0, . . . , k — 1. Kergin interpolating operators enjoy many interesting 
algebraic properties including the following. 

(1) /C[A] is independent of the ordering of the points in A, and 

(2) /C[B] o /C[A] = /C[B] for every B C A. 

In Theorem [231 Kergin operators are defined only for entire functions. 
Andersson and Passare [H [2], showed that Kergin operators JCd can ac- 
tually be defined on 0{D) where Z) is a C-convex domain in C^, i.e., the 
intersection of D with any complex line is connected and simply connected. 
In MJ^ this is simply ordinary convexity if we replace "complex line" by "real 
line." 

There are many results on the approximation of holomorphic functions 
by Kergin polynomials. We offer a brief sample. 

Let i^' C -D be compact and set /C„, := /C/)[A„] where, for n = 1, 2, 3, . . ., 
An = [Ano, . . . , Ann] C K. FoT D with C^— bouudary. Bloom and Calvi 
[TU] gave conditions on the array {A„}„=i 2,... so that /C„(/) converges to 
/ uniformly on ii" as n — )■ 00 for every function / holomorphic in some 
neighborhood of D. They utilized an integral representation formula for the 
remainder / — /Cn(/) proved by Andersson and Passare [1]. 

More in line with the ideas in this work, we call an array {A„}„=i 2,... 
extremal for a compact set K if /C„(/) converges to / uniformly on K for 
each / holomorphic in a neighborhood of K. For K C M"', Bloom and Calvi 
[TT] proved the following striking result. 
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Theorem 24. Let K C M'^, d > 2, be a compact, convex set with nonempty 
interior. Then K admits extremal arrays for Kergin interpolation if and 
only if d = 2 and K is the region bounded by an ellipse. 

Thus, for example, the real disk 

B = B2 = {{zi, Z2) E : Im^i = lmz2 = 0, (Rezif + (Re^s)^ < 1} 
admits extremal arrays for Kergin interpolation. 

11. Open problems. 
We conclude with some open questions, a subset of which comes from 

0- 

(1) Is the converse of Proposition [6] true? 

(2) Does (2) =^ (4)? 

(3) If an array lies in the Shilov boundary Sk of K, does (4) =^ (3)? 
The Chebysev-radii Bos array in B2 described in section 8 might 
give a counterexample. 

(4) Construct an explicit array in the ball B2 in section 8 satisfying (2), 
or, even better, (1). 

(5) Find an example of a compact set K C C^, d > 1, for which one 
can explicitly construct Fekete points. 

(6) Do multivariate Leja sequences satisfy (1)? (4)? 

(7) One can define multivariate weighted Leja sequences; starting with 
any point Xi G K, having chosen Xi, ...,Xm G K we choose Xm+i G K 
so that 

\W{Xi,...,Xm,Xrn+l)\ = max\W {Xi, Xm, 

Do these yield asymptotic weighted Fekete arrays? 

(8) For K C C''' compact and L— regular, does there exist c = c{K) > 1 
such that Fekete arrays of order en, n = 1,2,, form an admissible 
mesh (AM) for K7 

(9) For K G compact, L— regular, and polynomially convex, if a 
triangular array satisfies {Ga} is 9 — aT for K, is (4) satisfied? Is 
the converse true? Note if ci = 1 this equivalence is (essentially) 
Theorem [2l 

(10) Let K C C^, d > 1, be L— regular. If one takes asymptotic Fekete 
points, can the corresponding polynomials be used to recover the 
pluricomplex Green function Vk as in Theorem (TUP 
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